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We derive and discuss the temperature dependance of the condensate and noncondensate density 
profile of a Bose-Einstein condensate gas with Feshbach resonance in a parabolic trap. These quanti- 
ties are calculated self-consistently using the generalized Hartree-Fock-Bogoliubov (HFB) equations 
within the Bogoliubov approximation. At zero temperature, the HFB equation can be solved by 
means of a variation method that give the low excitation spectrum. Moreover, within the two-body 
collision theory, we estimate the relationship between the atom number and the external magnetic 
f"^ ■ field B, it is in good agreement with the data in recent experiments. 
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Bose-Einstein condensates [1-4] of atomic gases offer new opportunities for studying quantum-degenerate fluids. 
Almost all the essential properties of Bose-Einstein condensate systems such as the formation and shape of the 
condensate and noncondensate, the nature of its collective excitations and statistical fluctuations, and the formation 
and dynamics of vortices are determined by the strength of atom-atom interactions. In contrast with the situation in 
traditional superfluids, the strength of inter-particle interactions in the atomic condensate can vary over a wide range 
of values[5-12]. In particular, the scattering length that characterizes the atom-atom interactions can be negative, 
Q\ , corresponding to an effective inter-atom attraction. 

Most recently, in trapped atomic Bose-Einstein condensation, Ketterle's group reported evidences for modifying the 
scattering length by magnetic-field-induced Feshbach resonance[13]. Feshbach resonance have been studied twenty one 
years ago [14] at much higher energies, but the Feshbach resonance energy observed in ultracold atoms can be tuned 
to near zero. The theoretical studies of the ultracold atoms with Feshbach resonance [15, 16] show that the two-body 
interactions responsible for the Feshbach resonance produce an additional condensate of molecules [17], which differs 
qualitatively from the properties of a single condensate. In this letter, we discuss the temperature dependance of the 
hybrid atomic/molecular condensate and noncondensate density profile as well as the excitation spectrum. Especially, 
for an inhomogeneous system with negative scattering length, the excitation spectrum shows a upper constraint on the 
atom number, which, together with the conventional estimation of the condensate atom number, gives the relationship 
between the atom number and the external magnetic field B. The theoretical results are in good agreement with the 
experiment. 

The binary atom Feshbach resonances studied by varying a strong external magnetic field in an alkali-atom trap are 
hyperfinc-induccd spin-flip processes that bring the colliding atoms to a bound molecular state of different electron 
spin. This process can be described by the Hamiltonian 
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Ctf ; H FR = a / d 3 r$+ (r)¥ a (r)# B (r) + h.c, (1) 



where ^ m (r), ^+ (r) (\t a (r), $£(r)) are the annihilation and creation field operators of the molecules(atoms), a 
stands for the transition matrix element proportional to the overlap of the molecular continuum and bound state 
wave functions. Usually, the transition matrix element depends on the magnetic field as a ~ ^/A a A 2 /2\B — Bq\, 
where Bq is the resonant magnetic field and A characterizes the width of the resonance as a function of B. The 
Hamiltonian Hfr together with the atomic Hamiltonian 

H a = J d 3 r*+(r)[-^ + V a (r) - f* a ]* a (r) + d 3 r^(r)^(r)* Q (r)* Q (r), (2) 

the molecular Hamiltonian 

H m = f d 3 r*+(r)l-^M + ~ »™ + e ^™(r) + ^ f d?r^>+ (r)*+ (r)* m (r)* m (r), (3) 
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and the atom-molecule interaction Hamiltonian 



H am =\f d 3 r*+(r)*+ (r)* m (r)*„(r) (4) 

forms a total Hamiltonian H = H a + H m + H am + Hfr, which governs the dynamics of the system under investi- 
gation. Here, V a i m \{r) represents the trapped potential for atom (molecule), Ma(m) is the chemical potential of the 

atoms (molecules), \ a ( m ) = -757^ , with M being the atomic mass and a a ( m ) the s-wave scattering length of the 
atom-atom interaction, A denotes the coupling constant of the atom-molecule interaction, and e is the energy of the 
intermediate molecular state relative to the continuum of the incident atoms. 

Separating out the condensate part in the usual fashion (Bogoliubov approximation) i.e. 

*a(m)(0 = 4>a( m )(r) + <fia(m){r), (5) 

where (f> a ( m )( r ) — (^a(m){ r )} plays the role of a spatially varying macroscopic Bose field of the atoms(molecules). The 
possibility that the resonantly formed quasibound atom pairs form a molecular condensate was previously suggested 
by Timmcrmans et al.[lS\. Using a Raman photonassociation process, the quasibound pairs start to be formed from 
the atomic condensate and form a molecular condensate [17]. Here, we assume that there are a large number of 
atoms and molecules in the condensate. It is easy to show that the operators (p a (m){ r ) an d V^ m ){ r ) obey the Bose 
commutation relations 

W a{m) (r)^ + a(m) {r')]=5(r-r'). (6) 
In terms of (p a {m) an d 4> a (m)i the Hamiltonian can be expanded as 

h = h q + h\ 

H = J d 3 r{r a (-^ - Ma + V a (r))4> a + y «0a0a 
V 2 A 

+ (~ 4]^ - Mm + V m (r) + e)(f) m + -^0 TO m <Am0m 

+ Ha<Pa<P* m <Pm + OL(j)* m (j) a (j> a + a<Am0a0a} ( 7 ) 



/ f V 2 A 

H = j d M^a (-7TTT- - Lla + V a (r))ip a +2\ a ip+(f> a (f>* a (t> a + -j(lfi+ip+(j) a (f) a + ip a ip a (f>* a (j>* a 



V 2 A 

+ ^m(-^ + £ - Mm + V m (r))(p m + 2\ m ip+ip m (j)* m (f> m + -Y^ti^tn^m^m + <Pm<Pm<i>m < K, 
+ Wi<Pa<l>m<l> m + ty} n <Pm ( f>*a ( t>a + a((p a $ a (f>* m + (p+ <f) m )} . 



(8) 



In derivation of eq.(7,8), the following coupling equations are used[15], 

V 2 

V 2 

^~4M + Am '0 m ' 2 + V "ii r ) + 6 + A|0 a | 2 }0m + Ot(j) a <j) a = fi m (f> m . (9) 

This coupling equations may be yielded by the expectation value of the Heisenberg equations 

im a = [^ a ,H],im m = [* m ,H], (io) 

and replacing the time derivatives by the chemical potentials 

ih(j)a — > Va<f>a,ih<j>m ~ » Mm 0m- 

The chemical potential of the molecules is twice the chemical potential of the atoms, in accordance with the condition 
for chemical equilibrium. The a— terms that couple the equations describe tunneling of pairs of atoms between <f) m and 
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<p a fields, it leads to the form of a second condensate — molecular condensate in an atomic Bosc-Einstein condensate [15- 
17]. Using the coupling equations(9), Timmermans et al.[l5] investigate the behaviors of the hybrid atomic/molecular 
condensates near- and off- resonance. The Hamiltonian (8) can be diagonalized by using the Bogoliubov transformation 



<Pa(r) 
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-Vj(r)0j 



(11) 



where the qusiparticle operators aj,ctj,fij, [ilf obey boson commutation relations 

[ai, at] = Sij, [ai, aj] = [at, at] = 0, 

[A,/3t] = ^.,[ / 3 i ,/3 j ] = [/3t,/3t] = o, 

[a+ p+} = [a i ,p j } = [a i ,pf]=0, 
and Uj(r),vj(r), Xj(r),yj(r) are c-number functions. Substituting eq.(ll) into eq.(8), one yields 

H ' = E^s + ^+E e ^-E^/ ^i^wi 2 -E e » / ^mi 2 



(12) 



with 



V 2 



V 2 



V 2 

^~4M + 6 + A '^ a ' 2 + 2Am ^" 1 ' 2 ~~ Mm + ^W)^ ~ X ra<Kn ( KnV3 = e 3 X 3' 

V 2 

^~4M + 6 + A '^' 2 + 2Am ^" 1 ' 2 ~ ^ m + ^W)^ ~ X ™<t>* m (f>* m x j = - e jVj- (13) 

In order to study the temperature dependence of the excitation spectrum as well as the spatial distribution of the hybrid 
atom/molecular condensate and noncondensate, we need to solve the coupled mean-field Bogoliubov equations(13), 
and the condensate equation(9) self-consistently. The calculations procedure can be summarized for an arbitrary 
confining potential as follows: First of all, we solve eq.(9) self-consistently, once <f> a and <f> m are known, the solution 
of Uj,Vj,Xj and yj can be generated. To illustrate this procedure, we present its first step of calculations analytically. 
The trapped potential considered here is taken to be an isotropic harmonic potential V a ( m )(r) = ^Mu> 2 ^r , for 
which 4> a and <f> m are spherically symmetric functions, 

0a(m)M =floo(r)Foo(0,VO, (14) 

with 

R 00 (r) = a*l\[^exp[-\a 2 r 2 ],Y m {6,i>) = -L,a = ((2)M^) 1 / 2 . 



Rather than solving the coupled equations (13) directly, we introduce a new method based on the auxiliary functions 

U 3 = A 3( r \j)a> V 3 = B j{ r \j)a, 
X 3 = C 3( r \j)m,yj = Dj{r\j) m , 

(15) 
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where \j) a ( m ) is defined by 

V 2 1 
[~ 2 (4)M + V)( r )lb')«H = hu a (m)(j + -j)\j)a{ m )- 

The reason for such selection is that the level shifts caused by atom-atom interactions weakly depend on the shape of 
the wave function. A combination of Eqs. (13-15) gives 

(hu) a {j + i) + A|0 m | 2 + 2A a |^ a | 2 - nJAj - MIK + 2a<p* m )Bj = EjAj, 



1, 
2' 

1, 



(hu a (j + -) + X\^ m \ 2 + 2X a \^ a \ 2 - n a )Bj - M*J* a + 2a<f>* m )A 3 = -EjBj, 



{hu m (j + -) + e + \\4> a \ +2X m \(p m \ - fJ, m )Cj - \ m ^ m ^* m Dj = ejCj, 
{hcj m (j + i) + e + A|0 a | 2 + 2A m |0 m | 2 - n m )Dj - \ m 4>* m 4>* m Cj = -ejDj, (16) 
the eigenfunctions and the corresponding eigenvalues are given by 

i± _ r 1 i* D±r„\ _ „±t„ A\r>± 



7 ± (r,i) 

f L_ 



D f = L±,.^ 1 ]',gf(r)=g ± (r,j)i?f, (17) 



and 



£±(r) = ±{(A«CC + 2«C) - [^a(j + \) - Ma + A|^„| 2 + 2A Q |0 a | 2 ]}, 
ef(r) = ±{\ m ^ m 4>* m - [huj m (j + -) - p m + e + 2A m |0 m | 2 ]}. (18) 



Here, 



and 



g ± (r,j) = -Wmfe 

huJmU + \) - Mm + e + 2A m |0 m | 2 - (e^) 2 

These explicit solutions enable us to construct the one-body density matrix 

p(r, r ) = p a (r, r ) + 2p m (r, r ), 
p a (r,r) = <t>* a {r)<t>a{r ) 

OO 

+ E K(rK(r')F? + <*(r)<(r')(l + Ff)], 

p=±,i=l 

Pm(r,r) = <j)* m (r)<j) m {r ) 

OO 

+ E K*W(r')/f + wTW»f(r')(l + /?)]. (19) 

p— ±,i— 1 

where = exp ( /3 1 E P' ) _ 1 and /f = exp ^ e p pi arc the Bose distribution for the quasiparticle excitations with energies 

i?f and ef , respectively. Setting r = r , eq.(19) follows the resulting particle density. 

We need to point out that eqs. (17) and (18) are results of the first step of the numerical calculations. To complete 
numerical calculations, we should repeat the above procedures until the eigenvalues Ej and ej do not depend on 
position r. In what follows, we present a variation method to study the excitations at zero temperature. This method 
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was first introduced in Ref. [18] to study the BEC ground state in the harmonic trap of boson system, and it was 
generalized in Ref. [19] to investigate the excited states in BEC. Considering eq.(13) as well as 

J [uj{r)u*{r) - Vj{r)v*{r)]dr = 1, 

and 

J [xj(r)x*(r) - yj{r)y*{r)]dr = 1, 
which were derived from the Bose commutation relation (6) we arrive at 

Ej = J U * 3 ( r )(-J^ + MM 2 + 2X a \^ a \ 2 + V a {r)) Uj dr 
+ J v*(r)(-^ + M<f> m \ 2 + 2K\4>a\ 2 + V a (r))v 3 dr 
- J u*{r){\ a 4>* a 4>l + 2a4>*Jv j (r)dr - J v*{r)(\ a cp* a d?* a + 2acf>* m )u j (r)dr, 

f V 2 

ej= / x*(r)(-— + e + \\0 a \ 2 + 2\ m \cb a \ 2 + V m (r))x J dr 
+ J y*(r)(-^ + X\^ a \ 2 +e + 2\ m \^ m \ 2 + V m (r))y J dr 
-J x*{r)\ m (j)* m (f)* m yj{r)dr - J y*j{r)\ m 4>* m <f>* m Xj{r)dr. 



(20) 



(21) 







v*(r) ) 








V*ji.r) ) 





For simplicity, we study only the case of the spherical harmonic trap. In this case, we may choose the trial wave 
functions of the excitation components Uj(r),Vj(r),Xj(r) and yj(r) in the form of the spherical harmonic oscillator 
wave function £„ r ,z, m with quantum numbers (n r , I, m): 

(u nr im,r), (22) 

where cj nr im is an adjustable scaling factor of variation. Eqs (20) and (21) show that Ej and ej take a similar form, 
hence we here discuss branches Ej of the excitation spectra in detail. For (n r , I, m) = (0, 1, 0), we have 

ov"' 

The excitation spectrum in this case is reduced to 



E = E[v,lu 010 ] = (1 + 2v 2 )[-huj QW + -h^] 

4 4 WoiO 

2M 3/2 u ow 5/2 



+ A(l + 2, 2 )A^/ 2 [ s , 

+ [2A a (l + 2v 2 ) - 2X a vVTT^]N a ^J 2 [^ 2 [ "°"> f' 2 

- 4avVTT^[— ^f^NU 2 ^^', (23) 
^oio + w m Tin 



where a 2 0W = M^sua. Similarly, for (n r ,l,m) = (1,0,0), we have 



6,0,0 = "ioo[^] 1/2 (^ - aUr 2 )re-^/2 Yo o{e ^ 
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And 



7 7 lo 2 

E = E[v,u lm ] = (1 + 2v 2 )[-hw lm + -Jl- ' 



where 



' t J 

4 wioo 

+ A(l + 2v 2 )N m ^J l 2 [^-f 2 f{uj WQ ,2u Jm ) 
+ [2A a (l + 2v 2 ) - 2\ a v^l + ^]N a w 3 J 2 [^f/ 2 f{^ wo ,oj a ) 

AavVTT^fi^^N^Wj^r/^ (24) 

/(*. y) - §(-^) 3/2 - 3(^) 5/2 + §(-^) 7/2 . 

2 x + y x + y 2 x + y 

Minimizing the energies of eqs (23) and (24) with respect to the variation parameters v, woio and u>ioo, we can 
determine the excitation spectrum for the mode (0, 1, 0) and (1, 0, 0). The numerical results are illustrated in fig.l and 
fig. 2. The dashed lines in figures show the excitation spectrum in atom BEC, i.e., a = X = 0. In contrast, the solid 
line are those of hybrid atomic/molecular condensates near Feshbach resonance. From fig.l and fig. 2 we see that while 
the excitation frequency for modes (0,1,0) increases due to the Feshbach resonance effect, the excitation frequency 
for the mode (1,0,0) decreases. We would like to point out that the numerical results presented here depend on the 
coupling constant as well as the parameter a. In Fig.l and Fig. 2, we let a = 5A a , and A a = 0.1 (arbitrary units). The 
other parameters in Fig.l and Fig. 2 are co m = lAuj a = 7500Hz, N a = N m = 10 6 . 

For clarity, we illustrate the above somewhat formal discussion by considering the binary atom system for a uniform 
system (V^( m )(r) — > 0), in this case Hfr gives a resonant contribution to the atom-atom interaction strength a a : 
a e ff = a (l + g~rg ) i where a is the off-resonant scattering length, and A characterizes the width of the resonance. 
For small p, these excitations are phonons, and their energy tends to zero with p. Hence, 

lia = M^ m \ 2 + KM -2acj)* m , (25) 

which leads to 

E 2 = E 2 {p) = (^) 2 P 2 (P 2 + 16nna eff ). (26) 

For a uniform dilute Bose gas with negative scattering length a e //, eq.(26) implys an instability of those modes 
with p 2 < 167rn|a e //|. For a gas in a trap, however, the wavenumber cannot be arbitrarily small, and the minimum 
value is of order p m i n — n/Ro{ Rq is the mean size of the ground state). Hence the system can remain stable if 
jtz > 167rn|a e // 1. Since the density is of order n ~ N/R 2 , this means that the critical number of the system is 



16 \a ef f\ 



For a positive scatering length a e //, however, there are not any constraints in N. The Bogoliubov quasiparticle theory 

V ■ 



shows that the condensate atoms N depends on the scattering length and satisfies (for (a e //y) 3 << 1). 
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*o = ^1--^-^). (28) 

The numerical results of eqs. (27) and (28) are illustrated in Fig. 3, which shows the atom number vs. external 
magnetic field B. The parameters in Fig.3 are N/V = N/R 2 = 10 15 / 'cm 3 , A = O.OlmT. 

To sum up, we have derived a set of four coupled equations of the atomic and molecular excitations within standard 
Hartree-Fock-Bogoliubov approximation. As shown in eqs. (9) and (13), the a terms describing the process that 
converts atoms into molecules play an important role in atomic/molecular Bose-Einstein condensation. In particular, 
two low excitation spectrum have been given at zero temperature, which show that the interaction between the 
hybrid atomic/molecular BEC increase one excitation mode, while they decrease another excitation mode. The mode 
(0,1,0) comes from the density fluctuation of the condensate like vibrating oscillation, in this sense that the mode 
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(0,1,0) increases near the Feshbach resonance indicates the presence of the Feshbach resonance enhance the density 
fluctuation like vibrating oscillation in atomic/molecular condensation system, whereas the breath mode(like breathing 
oscillation) (1,0,0) decrease near the Feshbach resonance. Within the two-body collision regime, we show the atom 
number remained in BEC vs. the external magnetic field B, the result is in good agreement with the recent experiment. 
This work removes from consideration of the case at resonance, since at resonance the Bogoliubov approximation is not 
available( at resonance, there are few atoms in condensate). The contributions of the noncondensate atoms (molecules) 
to the excitation spectrum is also ignored (see eq.(9)). These need further investigations. 
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Figure Captions 

Fig.l: Excitation spectrum of Mode (0, 1,0) vs. number of atoms. Dotted and dashed line indicate those with and 

without Feshbach resonance, respectively. 

Fig.2:Same as fig.l. But for mode (1,0,0). 

Fig.3:The number of atoms in condensate vs. magnetic field B. 
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